Confinement effects in antiferromagnets 
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Phase equilibrium in confined Ising antiferromagnets was studied as a function of the coupling 
(v) and a magnetic field (h) at the surfaces, in the presence of an external field H . The ground 
state properties were calculated exactly for symmetric boundary conditions and nearest-neighbor 
interactions, and a full zero-temperature phase diagram in the plane v-h was obtained for films 
with symmetry-preserving surface orientations. The ground-state analysis was extended to the H-T 
plane using a cluster-variation free energy. The study of the finite-T properties (as a function of 
V and h) reveals the close interdependence between the surface and finite-size effects and, together 
with the ground-state phase diagram, provides an integral picture of the confinement in anisotropic 
antiferromagnets with surfaces that preserve the symmetry of the order parameter. 

PACS numbers: 75.70.-i; 68.35.Rh 
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I. INTRODUCTION 

Confinement effects play an important role in the ther- 
modynamics of several materials such as polymers, liq- 
uid crystals, and magnets. For example, capillary con- 
densation stands as a well known example of how phase 
equilibrium is affected by the confluence of surface and 
finite-size effects. In particular, due to the wall-particle 
interaction, a fluid between two plates undergoes a gas- 
liquid, ira^jsition at a lower pressure than it does in the 
bulk.EfHa'El These effects of confinement are due to the 
additional contributions to the thermodynamic potential 
of the solvation force (finite-size effect) and the wall-fluid 
intcrfacial tension (surface effect). eI 

A more complicated physical situation arises in the 
case of thin films of polymer mixtures on selective 
substrates.B An AB polymer mixture which undergoes 
a phase separation below a bulk critical temperature Tc, 
develops, when cast into a thin film over a substrate, 
an interface between the two phases which runs parallel 
to the substrate. This interface appears provided there 
is a substrate affinity for one of the components — the 
confinement is established between the polymer-air and 
polymer-substrate boundaries. 

A model fluid confined between two parallel walls that 
exert opposite surface fields, has been often considered 
in order to investigate the uadMr hr i inffl phy si cw hi . aTOtc ms 
with competing boundarics.cllal3'Blliyy'y yyizl In 
this case, the interplay between wetting and phase sep- 
aration is very important, unlike the case of capillary 
condensation in which wetting plays a small role. The 
competition between surface effects leads to an inter- 
esting and unusual behavior: Phase coexistence is re- 
stricted to temperatures below the wetting tempera- 
ture Tu, even in the limit of infinite separation be- 
tween the plates. The wetting temperature depends 
on the surface field jand it can be far from the bulk 
critical temperatureO The aforementioned scens 
first predicted using a mean-field approximation! 
has been jeoa&yrBejd subsequently via Monte Carlo 
simulationglJil3ll3'Eil and transfer-matrix calculations in 



two dimensionsJll However, when the effect of gravity is 
considered phase coexi^tjipjce is restored up to the bulk 
critical temperature. EjOEJ 

The confinement studies described above deal with 
phase separating systems, in which the phases coex- 
isting along a line of first-order transitions have the 
same symmetry, e.g., ferromagnetic thin films. Sur- 
face effects in systems with ordering (antiferromag- 
netic) interactions have been investigated mostly within 
the context of hiaajjiL^jaUoys undergoing a first-order 
phase transition,EHlElt30o with particular emphasis 
on the surface-iuduced order and surface-induced disor- 
der phenomenaJ£3'Lj although some investigations-.! 



been done in the context of multilayer adsorptionn!ln3'E£l 
More recently, attention has turned to the surface crit- 
ical behavior of jta naj} !. a i llQyS ] [displaying continuous or- 
dering reactionsBil'EHlEilOHEl'E^ and, in particular, to 
the dependenee of the universality class on surface 
orientation rWA 

In this paper we investigate the interplay between 
finite-size and surface effects in Ising antiferromagnets 
in the presence of an external field. In particular, we 
are interested in systems with surfaces that preserve the 
symmetry of the order parameter. In other words, we 
shall study thin films which develop antiferromagnetic 
(AFM) ordering in each plane parallel to the surfaces. 
Our layered system can be described by the following 
Hamiltonian: 



H = Jb ^ cTiaj + Js ^ (T 



ij Gsurf 



H ''^-ih + H) a,, (1) 



iebulk 



zGsurf 



where the spin variable takes the value of -1-1 or —1 
depending if the spin at site i is pointing up or down, 
respectively. We have assumed that surface sites, in lay- 
ers 1 and N for an iV-layer film, experience a surface 
field h in addition to the external magnetic field H. On 
physical grounds, it is natural to expect that the pair 
interactions at and near the surfaces differ from those 
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TABLE I: Energies for the different ground-states discussed in the text. Nomeclature is as foUows; structure T / Ti / T nieans 
that both surfaces are ferromagnetic (parallel to the applied field) and the remaining N — 2 inner layers are antiferromagnetic. 
Structure 5' is a special case in which the surface layers are ferromagnetic (t), the subsurface layers are also ferromagnetic 
but in the opposite direction (J.) and the rest (A'^ — 4) are antiferromagnetic (ti). Bulk and surface coordination number are 
denoted by z and z^. See the text for further explanations. 



Tag Structure Energy 

"1 i / i / i Hi = z, + {\z + H)[N - 2) + 2{H + h) 

2 n/i/n n2^-Zs + {\z + H){N ^2) 

3 T/i/T m = zs - Az^ + {\z + H){N - 2) - 2{H + h) 

4 i/n/i Hi^Zs-^ziN -2) + 2[H + h) 

5 Ti / Ti / Ti = -zs - \z{N - 2) 

5' T/i/Ti/i/T H„> ^ Zs + zo-2z^-\z{N -A)-2h 

6 T/Ti/T H6 = Zs~^z(N-2)~2{H + h) 

7 i/T/i H7 = Zs-4zi + {^z^H){N~2) + 2{H + h) 

8 Ti/T/n H8 = ~Zs + {^z-H)iN-2) 

9 T / T / T H9 = zs + {\z - H){N - 2) - 2{H + h) 



in the bulk. We approximate the position dependence 
of the pair couplings, by allowing the nearest-neighbor 
intralayer surface coupling (Jg) to differ from the bulk 
one (Jb). Here we restrict ourselves to case of Jh > 
(antiferromagnetic), but we allow to assume any real 
value. Also, we specialize ourselves in the case of local- 
ized symmetric surface fields, i.e., the field at each surface 
is the same and acts only at the surface sites. In the re- 
maining of the paper, the effective pair interactions, the 
surface field, and the external magnetic field (H) shall be 
expressed in terms of the bulk AFM coupling (Jf, > 0). 
The ratio of surface to bulk coupling is then denoted by v. 



Confinement effects in the order-disorder transitions 
for the particular case of v — 1 and h > have been 
reported previously.tJ In this paper, we give a full de- 
scription of the surface and finite-size effects in terms of 
the variables h and v. The ground-state properties of 
the Hamiltonian (^ are derived in Sec. l|. This zero- 
temperature analysis is used to identify the different se- 
quences of ground states displayed by the film as a func- 
tion of the external field H. Moreover, it is shown that 
for antiferromagnetic systems with symmetry-preserving 
surface orientations and nearest-neighbor interactions, a 
zero-temperature phase diagram can be drawn as a func- 
tion of V and h, for any value of the number of layers 
N and external field_i/. 



[I], we use a cluster 



In Sec. 

variation free energyEII to describe the finite-temperature 
behavior of the system as a function of surface variables 
V, h and the number of layers N. Particular attention is 
devoted to the analysis of the critical curve (in the H-T 
plane) for each one of the different regions of the zero- 
temperature phase diagram. We close with a summary 



II. GROUND-STATE PROPERTIES 

In the absence of surface and finite-size contributions, 
that is in the bulk, the Hamiltonian (pi) reduces to 



WbuU 



(2) 



of the important results (Sec. IV) 



For a two-sublattice antiferromagnet such as body- 
centered or simple cubic, the Hamiltonian (|^) has three 
different ground states as a function of the external field 
H: ferromagnetic (i) for H < —He, antiferromagnetic 
(ti) for —He < H < He and again ferromagnetic (T) for 
H > He- The critical field He, equal to the coordination 
number z [recall that all quantities in Eq. (|l|) as well as 
in Eq. (^ are normalized to Jb], determines the point 
where the critical curve Tc{H) meets the field axis. 

For the AFM thin films studied here [see Hamiltonian 
(|l|)], the possible ground-state (GS) structures are listed 
in Table | along with their corresponding energy. We 
considered only the case h > since the results for /i < 
can obtained straightforwardly from the symmetry prop- 
erties of Hamiltonian (|l|) . The nomenclature in Table B is 
as follows: structure number 4 corresponds to j / Ti / i, 
which means that both surfaces are ferromagnetic (i) 
and that the remaining {N — 2) inner layers are antifer- 
romagnetically ordered. Structure 5', a special case to be 
discussed later in the paper, has both surfaces in a fer- 
romagnetic state (T), the subsurface layers are ferromag- 
netic but with magnetization in the opposite direction 
(i), and the remaining [N — 4) layers are antiferromag- 
netic. 

We arrived at the set of GS in Table | as follows. Since 
only nearest-neighbor interactions are included in the 
Hamiltonian and the (uniform) surface field acts locally 
at surface sites, the presence of long-period superstruc- 
tures can be ruled out. A possible set of ground states 
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Magnetic field (H) 

FIG. 1: Schematic ground-state sequences VII (a), I (b), and II (c). The succession (a) to (c) represents the evolution in the 
ground-state sequence as we increase the surface pair interactions for neutral boundary conditions. The domain of stability of 
each structure is different from each other, but for h = the transition I ^ VII occurs at w = Vps = —zi/zq for which GS 5 in 
(b) has shrunk to zero width. The range of stability of GS 5 expands as we increase v. At v = Vm = {zo + z\)/z\ a transition 
between I in (b) and II in (c) occurs. As a reference, the value of the exchange energy for selected GS is indicated in the plots. 
See Fig. ^ for the evolution of sequences I and II with h. Also see the text for further explanations. 



for Hamiltonian (|lj) was then constructed by combining 
all possible surface and bulk ground states. For the sake 
of definiteness, let us consider a body-centered cubic film 
with surfaces in the (110) direction. The bulk ground 
states consist of two ferromagnetic structures (with op- 
posite magnetization) plus an ordered CsCl-type AFM 
structure. The (110) surfaces constitute face-centered 
rectangular lattices, for which the possible ground states 
are a checkerboard AFM structure and two ferromagnetic 
states of opposite magnetization. The nine ground-state 
structures obtained by combining the surface and bulk 
ground states are listed in Table |[ These structures are 
ground states of Hamiltonian (|l|) in the limit of weak cou- 
pling between the surface and the subsurface layers. For 
strong coupling between the surfaces and the bulk, we 
found onljuone additional ground-state structure — GS 5' 
in Table pa 

The Hamiltonian in Eq. (^ distinguishes between the 
pair interactions in the surface layers from the rest, thus 
allowing us to define the surface coordination number Zg 
as a function of the surface coupling parameter v 

Zs(v) = ZqV + Zl, (3) 

where the intralayer and interlayer coordination are de- 
noted by zq and zi, respectively. Recall that all quan- 
tities in Eq. (|l]) are given in terms of Jb and, therefore, 
Zb in (^) actually accounts for the surface energy. For a 
bcc(llO) film, zq — A and zi = 2, and the bulk coordina- 
tion number is z — zq + 2zi. 

Even in the absence of an applied surface field h, the 
surfaces are under the influence of a "missing neighbors" 
field hm{v) that arises from the disruption of the trans- 
lational symmetry perpendicular to the surfaces. This 
missing neighbors field produces an inhomogeneous mag- 
netization profile. Thus, with increasing external field, 
the surfaces may turn into a ferromagnetic state before 
the bulk does. Application of a surface field h = —hm{v) 
can restore the magnetization profile to the homogeneous 
condition. Note that the missing neighbors field depends 
on V, since hm is a measure of the difference between the 



environment at the surfaces and in the bulk [see Eq. (^]. 
The missing neighbors field can be written as 

hra = Zs- Z^ Z(^V - (zq -|- Zi). (4) 

We derive this value for the missing neighbors field later 
in paper, by considering the stability of the different GS 
structures as a function of v and h. 

A direct comparison between the energies 
Tii{v,h,H,N) for each structure i gives the ground 
state for every set of values of the thermodynamic 
variables (see Table ||). However, it is more useful 
and less tedious to consider a physical sequence of 
GS structures (as a function of the applied field) and 
examine its domain of stability as we vary the surface 
variables v and h. 

As a starting point, consider the following case: Upon 
the application of an external field H (in either direction), 
a film with v ^ 1 and h ^ will pass from an AFM state 
in all layers (small \H\) to a state with ferromagnetic 
surfaces and an AFM bulk and, finally, for large to a 
ferromagnetic state in all planes. This case is represented 
by the sequence 1-4-5-6-9 of GS structures [a schematic 
view is presented in Fig. |l|(b). See also Table | for the 
nomenclature]. The characteristic value of the external 
field at the transition between different GS structures is 
indicated in Fig. |l|. In general, the transition between GS 
structures A and B occurs at Hab, which is determined 
by equating the corresponding energies. 

Ground-state sequence 1-4-5-6-9 (hereafter referred as 
I) in Fig. [^(b), provides some useful insight on confine- 
ment versus finite-size effects. At the beginning of this 
Section we considered the ground states of an infinite 
antiferromagnet, which in the nomenclature of Table |, 
correspond to GS sequence 1-5-9 in the limit oi N ^ oo. 
Thus, ground-state structures 4 and 6 are due the con- 
finement effects. When either GS 4 or GS 6 become un- 
stable in favor of GS 5, surface effects are lost and the 
film is subject only to the finite-size effects. 

An external surface field will produce an asymmetry in 
the GS sequence since the Hamiltonian is not invariant 
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under the transformation ai — s- —ai, H — > —H. Applying 
a surface field h > reduces the surface ferromagnetism 
(I) in GS 4 and enhances it in GS 6 (t). The domain of 
stability of GS 4 shrinks to zero when Zg + h becomes z. 
This particular value of h defines the missing-neighbors 
field [Eq. d)]. 

Applying an external surface field is not the only way 
to eliminate surface effects in AFM thin films. A ho- 
mogeneous condition can also be attained in the film by 
setting neutral boundary conditions {h = 0) and increas- 
ing the pair interactions at the surfaces to a given value 
Vm- The characteristic value of the surface coupling that 
compensates for the missing neighbors effect is given by: 



{zq + zi)/zo. 



(5) 



For this value of the surface coupling GS 4 and GS 6 
become unstable simultaneously [Eq. (|5|) is equivalent to 
the condition z^ — z]. 

For values of the surface coupling larger than , keep- 
ing the neutrality at the boundaries, ordering becomes 
stronger at the surfaces than in the inner layers. This 
situation is represented in Fig. ^(c). It is worth noting 
that the GS sequence 1-2-5-8-9, hereafter referred as II, 
is stable not only in the case of /i = but for a range 
of values of u > v„i and h. We will return to this point 
later in the paper. 

Reducing the surface coupling make surface ordering 
less stable, until v reaches the characteristic value 



(6) 



for which GS 5 becomes unstable [vps in corresponds 
to the condition = 0] . The remaining ground-state se- 
quence 1-4-6-9 (hereafter VII) are depicted in Fig. ^(a). 
Sequence VII remains unaltered for v < Vps regardless 
the strength of v: The phase coexistence between spin- up 
and spin-down is regulated by the surface field h. Large, 
negative values of v increase the critical-point tempera- 
ture. In the alloy terminology, sequence VII represents 
the situation of a binary-alloy thin film with an ordered 
bulk coexisting with a surface miscibility gap. This will 
become apparent in Sec. Ill where we discuss the finite- 
temperature properties of Hamiltonian (|l|) . 

Sequences I, II, and VII (Fig. ^ were obtained by an- 
alyzing the stability of the corresponding GS sequences 
upon variations of the surface coupling v for neutral 
boundary conditions. As expected, a similar variation of 
GS sequences will appear as we increase the surface field. 
Consider for example sequence II in Fig. |l](c): Setting 
higher values for the surface field eventually overcome the 
ordering tendencies at the surfaces. Ground-state struc- 
ture 8 then becomes unstable and sequence II turns into 
the new 1-2-5-6-9 GS sequence (III) depicted inFig. |(a). 
The asymmetry of sequence III is interesting. For very 
negative values of H sequence III looks like sequence II 
(in the same range of 77), with long-range order dictated 
by the surfaces. On the other hand, for large positive 
values of H , sequence III looks like sequence I, for which 
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FIG. 2: Sequence II in Fig. ^(c) becomes sequence III in (a) as 
the surface field is increased. The transition II III occurs 
at /iii-iii = zqv — (zo + zi) when GS 8 becomes unstable (see 
Table O) . A further increase of the surface field h establishes 
1-2-3-5 -6-9 in (b) as the stable GS sequence (IV). Observe 
the appearance of GS 5' and the disordered gap between GS 
3 and GS 5'. The characteristic field between GS 2 and GS 3 
is _ff23 = Zs — 2z\ ~ h. See the text. 



the bulk is responsible for the AFM ordering. This sim- 
ilarity is due to the fact that sequence I evolves into III 
when the surface field increases beyond h ~ z — Zs(v) 
(missing neighbors field) for < u < Vm- 

The homogeneous antiferromagnetic thin film (GS 5), 
with constant energy for given v and iV, becomes rapidly 
unstable with increasing h. For sufficiently large /i, GS 5 
is replaced by another zero-magnetization structure, GS 
5' in Table ||, with energy given by 



Zs + zq- 2zi - \z{N - 4) - 2h. 



(7) 



For a given value of the number of layers N and the 
coordination at the surfaces, Ji^ is constant while 7^5' 
depends only on h. When the surface field reaches the 
value of 



tlll-IV 



zqv -I- (zo + ^i), 



(8) 



structures 5 and 5' have the same energy. A unique fea- 
ture of GS 5 and GS 5' is that they remain degenerate 
over a finite range of the external field H . From Eq. (H) 
and Fig. ^(a) we can see that the ground state for an 
AFM film at a surface field value given by (H) is a mix- 
ture of GS 5 and GS 5' for H S (— z, — zq)- On average, 
a scan in H will show a layer magnetization of -|-^ at 
the surfaces together with subsurface magnetization of 
— \ and AFM bulk (zero magnetization). Thermal exci- 
tations destroy this degeneracy between GS 5 and GS 5' 
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FIG. 3: Ground-state sequences for intense surface fields and 
negative v. Sequence V in (a) is obtained from IV in Fig. 
^b), when the surface AFM-phase in GS 2 becomes unsta- 
ble upon the reduction the pair interaction at surfaces. A 
coexistence between spin-up (GS 3) and spin-down (GS 1) 
magnetizations is established at the surfaces [cf. Fig. ^(a)]. 
For sufficient negative values of v and large h, a line of first- 
order transitions, ending at a critical point, occurs outside the 
antiferromagnetic critical curve. Sequence I in Fig. ^(c) turns 
into sequence VI in (b) in the same way as III becomes IV 
(Fig. H), that is, replacing GS 5 by GS 5'. A difference arises, 
however, since in this case GS 4 is adjacent to GS 5', and a 
coexistence line will appear mside the AFM region. See the 
text for details. 



in most of the interval {—z, —zq) in favor of GS 5, except 
near the ends, i.e. H ~ —z and H ^ —zq, where GS 5' is 
pinned by the onset of stability of GS 3 and the presence 
of GS 6. Traces of the ground-state degeneracy between 
GS 5 and GS 5' are observable at low temperatures. For 
the other structures listed in Table ^, a transition similar 
to 5 — > 5' does not occur, mainly due to the symmetry 
in the boundary conditions. 

Ground-state sequence III evolves into l-2-3-5'-6-9 se- 
quence (IV hereafter) sd h — /im-iv. The situation is 
shown schematically in Fig. ||(b) for h > /im-iv- Ob- 
serve that the appearance of GS 3 has established a dis- 
order gap between GS 2 (AFM surfaces) and GS 5' (AFM 
bulk). This behavior is unique in the sense that in all pre- 
vious cases the ordered domain was a compact interval 
in H. This characteristic brings some interesting fea- 
tures into the H-T phase diagram, such as the splitting 
of the filjn's critical Tc{H) curve into two distinct critical 
curves.E3 

Increasing the surface field does not change sequence 
IV into another GS sequence. However, for a large value 



of /i, antiferromagnetic order at the surfaces becomes un- 
stable upon reduction of the surface coupling, and IV 
changes into sequence V composed of GS sequence 1-3- 
5'-6-9 for w < [Fig. |(a)]. Observe in Fig. |(a) that 
GS 1 is now adjacent to GS 3. The difference between 
GS 1 and GS 3 resides at the surfaces, which have oppo- 
site magnetization. This situation is reminiscent to the 
one found in sequence VII [Fig. 0(a)] where the surface 
field regulates the surface phase coexistence between up- 
and down-magnetizations. In sequence V, however, the 
line of first-order transitions is located outside of the well 
defined AFM region composed by GS 5' and GS 6. 

Ground-state sequence V is stable for large h and neg- 
ative V. Previously, we found that VII is the stable GS 
sequence for v < Vps and low h. A transition between 
V and VII certainly occurs, although it is mediated by 
the GS sequence l-4-5'-6-9 (VI) [see Fig. |(b)]. Finally, 
sequences I and VI are separated by the GS 5 to GS 5' 
transition, which in this case occurs at lower values of h 
since v < 0. 

We have discussed the several GS sequences that ap- 
pear in confined antiferromagnets along particular paths, 
namely, we fixed h = and varied the surface cou- 
pling (Fig. m) or alternatively, we fixed v > Vm and 
increased /i (Fig. |^). In general, however, the transi- 
tion from one GS sequence to another does not occur at 
constant v or h. The relationship between the surface 
variables (w, h) at the different transitions between GS 
sequences (see Table defines the domain of stability 
of each sequence, from I to VII, in the plane v-h. The 
corresponding ground-state phase diagram is shown in 
Fig. IJ. The phase diagram is symmetric with respect to 
/i = 0, with the negative- /i region obtained by replacing 
spin-up with spin-down and H with —H in Fig. The 
zero-temperature phase diagram provides a good refer- 
ence frame to interpret some of the features reported in 
previous wopkr«jH-Jpkiary-alloy thin films with ordering 
interactions .E3E3e3c1I The ground-state phase diagram is 
also a valuable guide for the investigation of the finite- 
temperature properties of Hamiltonian to be carried 
out in the next Section. 



III. FINITE-TEMPERATURE PROPERTIES 

The finite-temperature properties of Hamiltonian (|l|) 
were calculated using the cluster-vadation method 
(CVM) in the pair approximation (PA).E!l For the two- 
sublattice antiferromagnets considered in this paper, the 
physical aspects of phase equilibri um .under confinement 
are well captured by the PA-CVMEaJlj For bcc(llO) films 
with neutral boundary conditions, a comparison between 
the PA and the tetrahedron approximation (TA) has 
shown that only, the quantitative aspects are improved 
with the latter .cil For a general exposition of the cluster- 
variation method, we refer the interested readrai.to..t]jie 
excellent reviews available in the literature J13c3'CJ'lSc3 

The order-disorder transitions are described in the 
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Surface coupling (v) 



FIG. 4: Ground-state phase diagram v-h for antiferromagnetic thin films. Table [[]^ contains the characteristic values of the 
surface field that describe the boundary between the different regions. Regions I-III define compact antiferromagnetic domains, 
while for region IV a disordered (ferromagnetic) gap intervenes between the AFM order at the surfaces and the ordered bulk. 
For large negative values of the surface coupling, a line first-order transitions ending in a critical point, appears in regions V-VII. 
The intra- and interlayer coordination numbers are denoted by zo and Zi , respectively. See the text for further explanations. 



usual manner by subdividing the bcc or sc lattice into 
two interpenetrating sublattices a and /3. The long-range 
order parameter in the fc-layer defined as 



(9) 



where 'ti^^^-j is the Q;(/3)-sublattice magnetization in the 
k layer. 



TABLE II: Characteristic values of h and v defined by the 
transitions between difi'erent ground-state sequences. These 
characteristic values define the domain of stability of struc- 
tures I to VII. See phase diagram in Fig. ^ and the text for 
additional details. 



Surface field/coupling 



hi.iii = -/in-iii = —zqv + {zQ + zi) 

/lI-VI = ^III-IV = ZQV + {zo + Zl) 

/iV-VI = Z() + Zl 

hyi-YU — 20 

Hl-II = (zo + Zi)/zo 

m-vn = -zi/zQ 



With reference to the GS phase diagram of Fig. ^, re- 
gions I-III display long-range order, either at the surfaces 
(GS 2 and 8) or in the bulk (GS 4^). With the exception 
of sequence I, the critical curvescll obtained in regions II 
and III show a distortion at high temperatures. Our re- 
sults for the critical curve in these regions, summarized in 
Fig. H, can be explained using the ground-state analysis 
discussed in Sec. l|. 

Phase diagrams in region I are virtually independent of 
the parameters v and h, as can be seen in Fig. g(a). This 
behavior can be attributed to the fact that the AFM or- 
dering in region I is primarily due to the inner layers [see 
Fig. 0(b) and Fig. ||. Thermal excitations can promote 
spin flip at the surfaces, resulting in a lower degree of 
ordering at surfaces relative to the bulk. In contrast, re- 
gion II is characterized by a strong AFM ordering at the 
surfaces coupled with the AFM bulk [see Fig. ^(c)] , thus 



preventing the formation of a (separate) surface critical 
curve. Instead, the H-T phase diagram shows an increase 
in the transition temperature and a broadening in the 
external field region for which the stable phase is antifer- 
romagnetic. A relative small asymmetry in the critical 
curve is observed, due to the fact that the surface field fa- 
vors the stability of GS 2 over GS 8 [Fig. |(a)]. Thus, the 
distortion in the phase diagrams associated with region 
II stems from the relative stability of two ground-state 
configurations with the same symmetry, i.e., GS 2 and 
GS 8. 

A higher asymmetry in the phase diagrams is expected 
in region III, since the critical-curve shape is dictated 
by the surface ordering for H ~ —(zs + h), and by an 
AFM bulk (with low surface ordering) for ~ z. The 
difference in symmetry of the AFM structures at each 
AFM:FM boundaries [see the GS sequence in Fig. ||(a)], 
allows the surfaces to drive the phase transition for fields 
close to the (negative) critical field value. One can see 
that the surfaces are developing their own critical curve, 
which unfolds as a 'shoulder' in the phase diagram for 
negative apphed field [see Fig. ||(b)]. Characteristics 
such as the maximum temperature of the shoulder or 
its extension in H, are controlled by the surface vari- 
ables V and h. The critical field between GS 1 and GS 
2 [Hi2 — —{z-u + h)] makes apparent that the exten- 
sion of the shoulder depends on the surface field. The 
maximum temperature in the shoulder is about wTsurf, 
where Tgurf is the Neel temperature of the corresponding 
surface antiferromagnet. Here, as in the rest of the pa- 
per, the relevant thermodynamic variables are expressed 
in units of the (positive) AFM couphng. Thus, in the 
PA a square lattice has a maximum critical temperature 
/cTsurf = 4/ln4 « 2.88. 

As pointed out previously, region IV is characterized 
by the formation of a disordered gap between two differ- 
ent ground states [see Fig. §(b)]. At finite temperatures 
and deep inside region IV, the surfaces develop their own 
critical curve well separated from the bulk antiferromag- 
netic region [see Fig. ra(a) showing the critical curves for 
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FIG. 5: (a) Typical critical curves for regions I (symbols) and 
II (solid line, v = 2.41 and h — 3.14). The various symbols 
correspond to different values of the surface variables {v, h) 
as shown in the inset. Note that in region I the shape of the 
critical curve is virtually independent of i; and h. (b) Phase 
diagram for region III (solid line), showing the development 
of a 'shoulder' as a signature of the incipient surface critical 
curve. The values for the surface variables are v = 1.5 and 
h — 11. The antiferromagnetic domain is a compact region. 
The phase diagram for a square lattice is shown as reference 
(dot-dashed line) and to illustrate the process of separation 
between the bulk and surface critical curves (c/. Fig. H). Both 
in (a) and (b) antiferromagnetic bcc(llO) films with A'' = 14 
were considered and solved in the pair approximation of CVM. 
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FIG. 6: For an intense field at the boundaries (h > /is), the 
surfaces decouple from the inner layers (bulk) and develop 
their own critical curve. In (a) the critical curves of 14-layer 
antiferromagnetic thin films are shown for ii = 1.5, /i = 14 
(circles) and v = 1.5, /i = 18 (triangles). The bulk critical 
curve showed no difference from h = 11 %o h = 18, hence 
only the former case is depicted. The solid lines represent, in 
the case of the surface critical curves, the phase diagram of 
a square AFM, appropriately shifted. The solid line in the 
bulk phase is the one associated to A'^ = 12 in region I. The 
splitting between the surfaces and bulk critical curves occur 
at Ta (temperature of splitting) and Ha (field of splitting), 
when the surface field reaches the value of ha . Part (b) shows 
a detail of the phase diagram of 100-layer AFM film at the 
very point of splitting. 



a 14-layer film with v = 1.5, /i = 14 (circles) and ft, = 18 
(triangles)]. Since the surfaces are weakly coupled with 
the bulk, the surface critical curve scales with v, i.e., the 
zero-temperature width of the AFM ordering is z^v and 
the maximum critical temperature is vkTsn-ci- 

Between the situation of unconnected ordered domains 
and the phase diagrams observed in region III, there is 
the case in which the zero-temperature disordered gap 
transforms, via thermal excitations, into a disordered re- 
gion in the H-T plane right inside the compact AFM 
domain. An increment in the surface field translates into 
an increment in the height of the disordered region. At 
h = hs the AFM region splits into the surface and the 
bulk critical curves [see Fig. ||(b)]. At finite tempera- 
tures, the splitting value of the surface field hg plays the 
role of /iiii-iv: for h < hg the ordered region is compact 
whereas for h > hg there are two unconnected critical 
curves. 

Expressing the free energy F in terms of the long-range 



order parameters (^, the conditions determining the lo- 
cus of the splitting point are given by: 



A = det 



drjk drjk' 



0, 



^=0 ^ = 
dT ' dh 



(10a) 



(10b) 



Equation( |10aD defines the critical temperature, at fixed 
external conditions (T, i7, h and N), when the second 
derivatives oLthe free energy are evaluated in the disor- 
dered state.E3 Since A < in the ordered state, at the 
splitting point {T — Ts) A is a concave function of the 
external field vanishing at the splitting value of the mag- 
netic field Hg. In a similar fashion, one can see that A is a 
convex function of temperature, becoming zero at T = 
[see Fig. ^(b)]. Thus, the splitting point is defined as a 
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FIG. 7: Splitting field //j, as a function of the surface coupling 
for AFM bcc(llO) films with = 14 and JV = 100. The 
minimum of Hs is due to the reentrance at low temperatures 
of the bulk critical curve. Inset: Splitting value of surface 
field ha as a function of the surface coupling u, for the case 
of N = 100 (circles). A least-squares fit (solid line) gives 
ha = 4.07t;-|-6.36. Compare this with /im-vi = 4:V+6 obtained 
in Sec. O for the boundary between regions III and IV. See 
the text for further details. 
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FIG. 8: A line of first-order transitions (thick line) appears 
in each of the regions V-VI for negative values of the sur- 
face coupling. In all cases, the coexistence is between up- and 
down-ferromagnetism at the surfaces. In all cases the bulk 
critical curve is not affected by the presence of the first-order 
transitions. Since the antiferromagnetic domain (shaded re- 
gion) is symmetric around H = only the left half is shown. 
The calculations were done in the PA-CVM for the following 
values of the surface variables: v = —1 and h = 2 (VII); h = 5 
(VI); ft = 9 (V). 



s add le point of A in the T and H variables. Conditions 
(lot) account for this. 

Using conditions (^) we determined the splitting value 
of the external field Hg as a function of the surface cou- 
pling for thin {N = 14) and thick {N = 100) films. The 
results are shown in Fig. |^ for the case of bcc(llO) films. 
The particular shape of Hs{v) can be understood as fol- 
lows: Since the height of the critical curve associated with 
the surfaces scales with v [see for example Figs. ||(b) and 
^(a)] and because of the reentrance of the bulk critical 
curve, for small v the point of contact (splitting) between 
the two critical curves is shifted to higher values of H. 
As we increase the surface coupling, the splitting point 
moves (clockwise) along the bulk critical curve, reach- 
ing a minimum in H and increasing again towards the 
saturation value. 

We found that within the PA the minimmn in Hs{v) 
is not very sensitive to the total number of layers. For 
bcc(llO), 77™™ occurs at w - 1.74 while for sc(lOO) the 
Hs is minimum at v ~ 1.2. Again, this can be explained 
by considering the different Neel temperature values for 
sc and bcc lattices. The ratio between the latter and the 
former is ~1.4 (PA), which is comparable to the ratio 
of the corresponding if™'" (~1.45). The behavior of the 
other quantities of interest can be inferred from Fig. |^. 
The most interesting part, however, is contained in the 
inset of Fig. ^, which shows hs as a function of the surface 
coupling V. A least-square fit gives h^'^'^ = i.Olv + 6.36 
which is almost parallel (and very close) to /im-iv in Eq- 
(^. For sc(lOO) films similar results were obtained and 
a linear fit for hs gives hf = 4.04w -|- 5.20. Thus, the 
process of splitting occurs within a narrow interval of h. 

Due to the equivalence between the Neel point and the 
critical point of a ferromagnet in zero field, the finite- 
temperature behavior of AFM thin films, as a function 



of the surface coupling v and H ^ h — 0, is equivalent to 
the multicritical phenomeaa occurring at the surface of 
semiinfinite ferromagnets.EJ In our case, negative surface 
pair interactions give rise to a line of first-order transi- 
tions in regions V-VII (see Fig. ^). In all cases the coex- 
istence line separates surface ferromagnetic phases with 
opposite magnetization, that have the same symmetry. 
The bulk, however, may have different symmetry at each 
side of the coexistence line, thus modifying the shape of 
the first-order line at finite temperatures. This can be 
observed in Fig. |^, where the surface coexistence curve 
is drawn w = -1 and ft, = 2 (VII), /i = 5 (VI), h = 9 
(V) . In each case, the coexistence curve ends in a critical 
point which is close, as expected, to the Curie point as- 
sociated with the (2D) surface lattice, i.e. ^ \v\Tc- In all 
the three regions V-VII, the AFM bulk remains undis- 
turbed by the presence of the surface coexistence line. 
At V = vn and H ^ h ^ the critical end point reaches 
the second-order critical curve at the Neel temperature 
T = Tjv. The multicritical behavior is the (trivial) su- 
perposition of two independent e«tical behaviors which 
do not interfere with each other .c2l 



IV. SUMMARY AND CONCLUSIONS 

In this paper we performed an analysis of the con- 
finement effects on antiferromagnets with symmetry- 
preserving surface orientations. The ground-state prop- 
erties of the model, an Ising Hamiltonian with nn-pair 
interactions in the presence of external bulk and sur- 
face fields, shows an interesting structure. A zero- 
temperature phase diagram in the surface variables v 
(surface coupling) and h (surface field) was obtained for 
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two-sublattice antiferromagnets. In this case there are 
seven different regions in the ground-state phase diagram. 
Each region is characterized by a particular sequence of 
ground states as a function of the external field. An 
analysis of the ground-state phase diagram explains (and 
sometimes even anticipates) some of the features found in 
the H-T critical curves. Together with an examination of 
the finite-temperature behavior in each of the aforemen- 
tioned regions, our analysis showed that the interplay be- 
tween the surface variables v and h defines the thermody- 
namics of confinement in ordering systems. For example, 
the splitting of the critical curve into surface and bulk 
contributions results from the simultaneous application 



of seemingly competing contributions v > Vm (ordering) 
and h > ftn-iv- At the other extreme, the development of 
a surface coexistence line for v < and h > represents 
a particular case of magnetic surface reconstruction. 
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